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The spectrum of meson and diquark excitations of dense cold quark matter is investigated in the 
framework of a Nambu-Jona-Lasinio type model for light quarks of two flavors. It was found out 
that a first order phase transition occurs when the chemical potential /i exceeds the critical value 
^ . jj,c = 350 MeV. Above fi c , the diquark condensate (qq) forms, breaking the color symmetry of strong 

' interaction. The masses of ir- and cr-mesons are shown to grow with the chemical potential fi in 

| the color-superconducting phase, but the mesons themselves become almost stable particles due to 

the Mott effect. Moreover, we have found in this phase an abnormal number of three, instead of 
five Nambu-Goldstone bosons, together with a color doublet of light stable diquark modes and a 
color-singlet heavy diquark resonance with the mass ~ 1100 MeV. In the color symmetric phase, i. e. 
for /j, < (j, c , a mass splitting of diquarks and antidiquarks is shown to arise if fi 7^ 0, contrary to the 
case of vanishing chemical potential, where the masses of antidiquarks and diquarks are degenerate 
at the value ~ 700 MeV. 
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I. INTRODUCTION 



One of the challenging problems of elementary particle physics is the investigation of hot and/or dense strongly 
(-h , interacting matter. At normal conditions (low temperatures and baryon densities), it is in the hardronic phase, where 
O ,' quarks and gluons are confined and chiral symmetry is broken. It is widely believed that at high temperature strongly 
interacting matter exists as a quark-gluon plasma (QGP). Another example of matter under extreme conditions is the 
interior of compact stars (low temperature and rather high baryon densities) which presumably consists of nothing 
else than dense cold quark matter. Hopefully, the properties of quark matter in extreme conditions become observable 
in relativistic heavy- ion collision experiments and/or through modifications of star evolution processes. 

Clearly, the underlying theory of strongly interacting matter, both in the vacuum and in extreme conditions, 
is QCD. Unfortunately, a proper and reliable quantitative description of quark matter in terms of a perturbative 
expansion of QCD is available only for asymptotically high values of temperature and/or chemical potential (densities). 
Nevertheless, perturbative QCD calculations show that at high temperature and low baryon density there is actually 
a QGP phase, with quarks and gluons being free particles and with the chiral symmetry being restored. Recently, 
a progress has been done in extending the lattice QCD approach to small nonvanishing values (< 140 MeV) of 
the chemical potential /1 (see, e. g. |l| and references therein). However, for the range of interest, i. e. for /1 ~ 
300 — 400 MeV, lattice QCD does not help. Anyway, perturbative QCD considerations, performed at low temperature 
and asymptotically high baryon density (large /S) , indicate the occurrence of a new color superconducting phase of cold 
quark matter 0- The confinement is also absent in it, but the ground state, unlike the case of QGP, is characterized 
by a nonvanishing diquark condensate (qq) . 

At moderate baryon densities, the QCD coupling constant is too large so that perturbation theory fails in this case. 
Obviously, the investigation of quark matter can then be suitably done in the framework of an effective quark model, 
e. g. in a chiral quark model of the Nambu-Jona-Lasinio (NJL) type Q including various channels of four-fermion 
interactions. Recently, on the basis of NJL type models, it was shown that a color superconducting phase might 
yet be present at rather small values of fi ~ 350 MeV, i. e. for baryon densities just several times larger than the 
density of the ordinary nuclear matter (see, e. g., 0,0 for a review). But this density is presumably reached in 
the cores of compact stars. On the other hand, the color superconducting quark matter inside compact stars, if it 
exists, can reveal itself through modifications of the star evolution process. The latter is the subject of astrophysical 
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observations and might be actually picked out from the data which are now being collected. By these reasons, the 
color-superconducting quark matter surely deserves a more detailed study. 

On the microscopic level, the processes running inside compact stars and in the fireballs created in heavy-ion 
collisions are understood as governed, to much extent, by strongly interacting quarks and gluons. Recent experiments 
at RHIC have shown that the hot and dense quark matter reached at the collider is far from full asymptotic freedom; 
instead, it is 'strongly coupled' @- In this case, the correlations between quarks and antiquarks, in particular due 
to composite mesons, are not negligible. This is why we start an investigation of the pion and cr-meson in the quark 
matter in extreme conditions. These particles are expected to be numerously produced in heavy-ion collisions. We 
also investigate diquarks, because they are important in determining baryon properties. 

In the present paper we work with an extended two-flavor NJL model in order to study the ground state of cold 
quark matter and some lightest meson and diquark excitations in the two-flavor color superconducting phase (2SC). 
For simplicity only a single quark chemical potential fi (common for all quarks) is used in the model. Moreover, we 
restrict ourselves to the region fj, < 400 MeV (notice that at higher values of /j, the color-flavor-locked phase is a more 
preferable one 3). Insofar as the phase diagram of quark matter has been discussed in a lot of other papers, we focus 
on the meson and diquark excitations only. Partially, we have already addressed this problem in our previous paper 
Q, where we had shown that in the NJL model under consideration an anomalous number of Nambu-Goldstone 
bosons is present when the initial SU(3) C symmetry of the model is spontaneously broken down to SU(2) C . Now, 
we are interested in the investigation of the masses of the remaining heavy diquark as well as of the pion and the 
cr-meson. 

The paper is organized as follows. In Section II, we present the Lagrangian of the extended NJL model as well 
as an equivalent Lagrangian that contains meson and diquark fields coupled with quarks. Using the Nambu-Gorkov 
formalism, we then derive an expression for the quark propagator, which is suitable for the study of the 2SC phase. 
Its poles provide us with the (anti)quark dispersion relations. Moreover, in this section, the gap equations for both 
the chiral and diquark condensates are derived. The meson and diquark masses are calculated in Section III; they are 
found to be monotonously increasing functions of fi in the 2SC phase. Moreover, the pion and cr-meson become almost 
stable particles in the 2SC phase due to the Mott effect (they are allowed to decay through electro- weak channels only). 
Finally, we discuss the mass splitting between the diquark and antidiquark states arising in the color-symmetric phase 
and the phenomenon of an anomalous number of NG-bosons in the 2SC phase. Section IV contains our conclusions 
and discussion. Some technical details are worked out in two Appendices. 



II. TWO-FLAVOR NJL MODEL 



A. Lagrangian 



In the original version of the NJL model 3] , the four-fermion interaction of a proton (p) and neutron (n) doublet 
was considered, and the principle of the dynamical breaking of chiral symmetry was demonstrated. Later, the (p,n)- 
doublet was replaced by a doublet of colored up (it) and down (d) quarks (or even more generally, by a flavor triplet) 
in order to describe phenomenologically the physics of light mesons 

II H M G3, diquarks [H [jj as well as the 
meson-baryon interaction 0,^1- In this sense, the NJL model may be thought of as an effective theory for low-energy 
QCD. 1 (Of course, one should keep in mind that, unlike QCD, quarks are not confined in the NJL model.) At present 
time, the phenomenon of dynamical (chiral) symmetry breaking is one of the cornerstones of modern particle physics. 
It has been studied, for example, in the framework of NJL-type models with external magnetic fields [l6| . in curved 
space-times [T^ |. in spaces with nontrivial topology [Tsl ]. etc. In par ticular, the properties of normal hot and/or dense 
quark matter were also considered within such models [Tol IT^. l20l l2l | . NJL-type models still remain a simple but 
useful instrument for the exploration of color superconducting quark matter at moderate densities 0, l22l |23| , where 
analytical and/or lattice computations in QCD are hindered. 

Instead of formulating a thermal theory in Euclidean metric (which is natural in statistical physics when one wants 
to derive a grand potential), we extend the Lagrangian for the two- flavor NJL model in Minkowski metric by the 
inclusion of the chemical potential fi and obtain 



L„ 



{qq) 2 + {qij 5 Tqf 



Go 



£ [■ 

A=2,5,7 



q c ij 5 T 2 X A q 



qij T 2 \Aq 



(1) 



Indeed, let us consider two-flavor SU(3) c -symmetric QCD. By integrating over gluons in the generating functional of QCD and further 
"approximating" the nonperturbative gluon propagator by a (5-function, one arrives at an effective local chiral four-quark interaction of 
the NJL type, describing low-energy hadron physics. Moreover, after the Fierz transformation of the interaction terms, one obtains an 
NJL-type Lagrangian that describes the interaction of quarks in the scalar and pseudo-scalar (qq) as well as in the scalar diquark (qq) 
channels (see, e. g., Lagrangian in this section). 
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where the quark field q is a flavor doublet and a color triplet as well as a four-component Dirac spinor, q c = Cg*, 
q c = q l C are charge-conjugated spinors, and C — ij 2 ^ is the charge conjugation matrix (t denotes the transposition 
operation). Here, the isotopic symmetry of quarks is implied (?71q = TOq = too) and the quark chemical potential 
fi > is the same for all flavors. Pauli matrices r a (a = 1,2,3) act on the flavor indices of quark fields, while 
the (antisymmetric) Gcll-Mann matrices Xa contract with the color ones; hereafter, the flavor and color indices are 
omitted for simplicity. Clearly, the Lagrangian L q is invariant under transformations by the color SU(3) C as well as 
by the baryon U(1)b groups. In addition, when the current quark masses vanish (too = 0), this Lagrangian resumes 
the (chiral) SU(2)l x SU(2)r symmetry (chiral transformations affect flavor indices only). 2 Finally, we remark that 
Lagrangian is C-even in the vacuum, i. e. it is invariant under the charge conjugation if fi = (q — > q c = Cq 1 , 
q — > q c = q t C), which is not the case for dense quark matter where the violation of C-parity is induced by a 
nonvanishing baryonic chemical potential. 

Throughout all our calculations, we assume that the model parameters, i. e. the ultraviolet cutoff A 3 , the current 
quark mass, and the coupling constants, do not change with /i. Their values are fixed in the vacuum so that the model 
can reproduce the experimental value of the pion mass Mj = 140 MeV, the pion weak-decay constant F v — 92.4 MeV, 
and the value of the chiral quark condensate (qq) = —(245 MeV) 3 . It has been shown in previous papers that a 
convenient set of parameters for the vacuum case is G\ — 5.86 GeV~ 2 , A = 618 MeV, and mo = 5.67 MeV, leading 
to the constituent quark mass 350 MeV. (One can follow, e. g., the parameter-fixing procedure explained in d,|2^| to 
obtain values close to these.) However, the definition of the constant G 2 that describes the interaction of quarks in 
the diquark channel is not quite transparent. This constant is not bounded by some experimental reason, and often, 
its value is fixed through a constraint (obtained after the Fierz transformation) that connects G 2 with constants in 
the quark-antiquark channels. Starting from the four-quark vertices provided by the one-gluon exchange, one obtains 
the constraint G2 — 3Gi/4 which we shall use in our further calculations. 

In principle, all constants in the diquark channels could be fixed if all interaction constants in the quark-antiquark 
channels were known. Unfortunately, only few of them are available from experiment, thereby preventing us from 
a unique definition of the constants in the diquark channel at the moment. There is still a freedom in fixing their 
values. This was the reason to use some model assumptions, such as the constraint that follows from the one-gluon 
exchange contribution. 



B. Quark propagator in the case of diquark condensation. Gap equations 



It is convenient to consider a linearized version of Lagrangian JQ) containing auxiliary bosonic fields which is given 
by the following form 



L 



l v id v 



M7° - o" 
iA* 



m - IJ TT a T a 



1 

4G7 



00 



4G 2 



-A*, A 



A [q C il b T 2 \Aq} 



■[qij 5 T 2 X A q C 



(2) 



As usual, the summation over repeated indices a — 1,2,3 and A, A' = 2, 5, 7 is implied throughout all our calculations. 
The equations of motion for the bosonic fields are 



a(x) = -2Gi{qq), 
ir a (x) = -2Gi(qij 5 T a q), 



A A (x) = 2iG 2 {q c i^T 2 X A q), 
A* A (x) = ~2iG 2 {qi 1 \ 2 \ A q C ). 



(3) 



Substituting into ©, one can easily obtain the initial Lagrangian Q). Just in this sense the two theories, the first 
one with L q and the second one with L, are equivalent. It follows from Q that the meson fields <r, 7r a are real, i. e. 
(<r(;r))t = <r(x), (n a (x))^ ~ n a {x) (the symbol f stands for the hermitian conjugation), whereas all diquark fields 
A^ are complex, (Aa{x))^ ~ A* A {x). Each A^ is an isoscalar (SU(2)l X SU(2)R-singlet). Moreover, all diquarks are 
Lorentz scalars and form an antitriplet (3 C ) fundamental representation of the color SU(3) C group, while the real scalar 
a and pseudoscalar ix a fields are color singlets. A nonvanishing value of the scalar diquark condensate, associated 
to a nonzero ground-state expectation value of some diquark field, (A^) ^ 0, breaks SU(3) C spontaneously down to 
SU(2) C (however, it does not violate the chiral SU(2)l x SU(2)r symmetry), whereas a nonzero expectation value of 
(a) ^ at too = indicates that the chiral symmetry is spontaneously broken. We assume hereafter that P-parity is 
conserved, i. e. (-K a (x)) = 0. 



2 Since Q = I3 + B/2, where ^3 = T3/2 is the third component of the isospin, and Q and B are electric and baryon charges, respectively, 
the Lagrangian is invariant under U(1)q transformations generated by electric charge, as well. 

3 There are divergent integrals in the model, and a regularization is needed. 

4 Some c onstants in the diquark channel can, however, be extracted from the nucleon mass, e. g. within the Bethe-Salpeter approach 

0m 
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Using the Nambu-Gorkov formalism, we put the quark fields and their charge conjugates together into a bispinor 
^ = (qa\ and Lagrangian (J2J, thereafter, simplifies to 



a ^-w?;.j i ... ( 4) 



4Gi 4G 2 2 V }C *' V 
where the following notations are adopted 5 

V+ = iYd v - m + - E, £>" = - m - M7° - £*, 

S = cr + i 7 5 7rf, S* = a + ij 5 TTT\ JC* = -A* a X a j 5 t 2 , IC = A a X a j 5 t 2 . (5) 

Now, let (er) = m — mo 7^ and (A a) 7^ for some A. Without loss of generality, one can always take advantage 
of the color symmetry of strong interactions and rotate the basis of diquark fields so that (A2) = A, (A£) = A*, 
(A5) = 0, (A7) = 0. As was shown in previous investigations, A = in the vacuum and for small /1, and the quark 
matter is thereby color symmetric. Once A acquires a nontrivial value, the SU(3) C symmetry is spontaneously broken 
down to SU(2) C and the 2SC phase is formed. According to the mean-field approximation approach, which we use 
here, the mean value of A should be subtracted from the diquark-quark vertices in (0J (or (J2J). Shifting the fields as 
follows: cr(x) — > a(x) + (a), A'^x) — » A^x) + A*, A 2 (x) — ► A2(x) + A, we absorb the mean values of a and A (with 
the exception of the quadratic terms) in the inverse quark propagator S 1-1 and obtain 

where er, 7r a , A^ and A A stand for the fluctuations around the mean values of mesons and diquarks rather than the 
original fields 6 , and 

c-i _ f - m + ^7°, AA 2 7 5 r 2 \ v _ f -a - ij S Tfr, A a j 5 t 2 X a \ f7 s 

\ -A*A 27 5 r 2 , vfdv - m - ^7° J ' \~A* Alb T 2 X Al -a-i^f 1 )■ {i > 

Now we can perform a series expansion, treating the term V, which contains only fluctuations of meson and diquark 
fields, as a perturbation. In the rest of our paper, we shall keep in mind the Feynman diagram rules for the calculation 
of two-point field correlators (Green's functions) in the 2SC phase, however, without drawing the corresponding 
graphs. Indeed, the term S^ 1 supplies us with the quark propagator S in the presence of a diquark condensate (qq), 
and the term V is responsible for quark-meson and quark-diquark vertices. In momentum space, the quark propagator 
( S S \ 

&(p) = I a ' a I (here p is the four-momentum of quarks) is represented by a 2 x 2 matrix, with respect to the 
\ J21, J22 J 

Nambu-Gorkov indices (in addition, it is a 2 x 2-, 3 x 3-, as well as 4 x 4 matrix in the flavor-, color-, and spinor 
spaces, correspondingly). After some algebra, one obtains for S i3 {p) ( 0, 0): 

S ij {p) = Slf{p)F Ig + S^ j (p)P h , (8) 

where P rg = diag(l,l,0), Pb = diag(0,0, 1) are diagonal matrices projecting onto the red-green and blue quark 
components in color space, respectively, and 

rg po + E- ~ po - E+ ~ I A+ A_ \ 

Sli = ^^7oA_ + ^7^7oA +! (9) 

ob _ 7qA- 7qA+ b _ 7qA- . 7qA+ qb _ Q b _ n rifV , 

Po - E p Q + E+ po - E+ Pq + E 




5 In the derivation of J1J, we used the well-known relations: c9* = —d v , Cy"C 1 = —(7")*, C"( 5 C 1 = (7 5 )* = 7 5 , t 2 tt 2 = — (t)*, 

, -i \ 

1 , ) ' 

6 One must not be confused with using the same notations both for the original fields and their fluctuations, because the original fields 
will never appear in the rest of the paper. 
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We used in eqs. an d (fTTTf) the projection operators 



A ± (P) = i (l ± 1^±^ , a ±(p1 = \ (l ± 

to separate the 'positive-energy' and 'negative-energy' parts of the quark propagator. Note the following useful 
properties of these projectors: 

7 A±7o = A T , 75A±7 5 = A±. 
Finally, the following notations have been used in eqs. 10 and (|1(J|) 



D + (p )=p 2 -(E+) 2 , D_(p )=pl-(E- A )\ E± = V(£±) 2 + |AP, E ± =E±^ (12) 



where E = yjp 2 + m 2 is the dispersion law for free quarks. The poles of the matrix elements and (I10J) of the quark 
propagator give the dispersion laws for quarks in the medium. Thus, we have E^ for the energy of red/green quarks 
and E~^ for the energy of red/green antiquarks. Below, we shall show that /i > m in the 2SC phase (see Fig. 1), and 
E can reach the value of [i. In this case, in order to create a red/green quark in the 2SC phase, a minimal amount of 
energy (the gap) equal to | A| at the Fermi level (E — /i) is required. Similarly, the energy of a blue quark (antiquark) 
is E~ (E + ), hence E~ = at E = //, and there is no energy cost to create a blue quark, i. e. blue quarks are gapless 
in the 2SC phase. (In contrast, to create an antiquark in the 2SC phase, we need the energy y/ (m + [i) 2 + |A| 2 for 
red/green antiquarks and the energy (m + (i) for blue one.) Given the explicit expression for the quark propagator 
S(p), we then calculate two-point correlators of meson and diquark fluctuations over the ground state in the one-loop 
(mean-field) approximation. 

The two additional quantities m and A are not free parameters of Lagrangian © . The constituent quark mass m 
is an indicator of the chiral symmetry breaking (if mo = 0, m vanishes when the chiral symmetry is restored). The 
gap A is related to the color symmetry breaking/restoration in a similar way. Both of them arc to be found from 
the requirement that the ground state expectation values of quantum fluctuations must be zeros. It is easy to see 
from © that in the mean-field (one-loop) approximation the two conditions (a) = and (A2) = are realized if the 
following two equations (gap equations) are true 

Here Ta* stands for the q-q-A^ vertex in the Nambu-Gorkov representation: Ta* = ( ^ \ ~ ) . In (|13|l and in 



-75T2A2, 

all similar formulas below, the calculation of the trace includes also a sum over Nambu-Gorkov indices, in addition to 
the spinor, color, and flavor ones. (Using the Feynman diagram terminology, one can say that the first term in each of 
eqs. Ijl3(l is a tree term whereas the second one represents the one-loop contribution.) After some trace calculations, 
we have from (|T3)l 

m-mp _ f d*q f E+ E~ 2E+ 2E~ \ 

2Gi " (2n)iE\ q 2 -(E+y + q 2 -(E-)2 + D + (q ) + D_(q Q )r ( j 

A f d^q ( 1 11 

4*A / -3-A — — + — — . (15) 



4G 2 J (2tt) 4 \D+{q ) D_(q ) 

To proceed further, wc employ the imaginary-time formalism, used in theories with finite temperature and chemical 
potential in Euclidean metric to obtain Green's functions. In these theories, at some finite temperature T, the 
integration over the energy- variable in each loop is replaced by a sum over Matsubara frequencies. The case of cold 
quark matter can be considered simply as the limit T — > 0. One should note here that we are interested in Green's 
functions in Minkowski metric rather than in Euclidean, and a continuation from one metric to the other is needed. 
Let us consider Green's functions as functions of energy only, forgetting, for a moment, about the three-momentum. 
At this point, we assume that after the limit T — > is taken, the resulting Green's functions are defined in a complex 
plane, and we associate the points lying on the imaginary axis with Euclidean metric. We continue the Green's 
functions to the real axis and consider the thus obtained new functions as being defined in Minkowski metric. Clearly, 
one should impose the constraint that such a continuation must reproduce the result obtained from quantum field 
theory in Minkowski metric when we put /i = 0. 

To calculate the integrals in (|14fl and i|15fl . we replace the integration over q by a sum over Matsubara frequencies, 
Lo n = (2n + 1)T, n = 0, ±1, ±2, .., followed by the limit T — > 0: 

9 :f(qo,q) — »iJjm,rV / -^4^f(i^ n ,q)- (16) 



- T^0 (27T) 
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FIG. 1: The constituent quark mass m (solid line) and the color gap A (dashed line) as functions of the chemical potential fi. 

Applying rule (|16|l to cqs. (|14(l and (|15|l . we obtain the gap equations for to and A in the cold quark matter (T = 0) 
(for more explanations, see Appendix A): 



Since the integrals in the right hand sides of these equations are ultraviolet divergent, we regularize them and the 
other divergent integrals below by implementing a three-dimensional cutoff A. 

The system of eqs. (|17H and (|18f) has two different solutions. As we have already discussed, the first one (with 
A = 0) corresponds to the SU(3) c -symmetric phase of the model (normal phase), the second one (with A ^ 0) to 
the 2SC phase. As usual, solutions of these equations give local extrema of the thermodynamic potential Jl(m, A) 
7 , so one should also check which of them corresponds to the absolute minimum of il. Having found the solution 
corresponding to the stable state of quark matter (the absolute minimum of Q), we obtained the behavior of the gaps 
to and A vs. chemical potential (see Fig. 1). The region /i < ji c = 350 MeV is the domain of color symmetric quark 
matter because il in this case is minimized by to 7^ and A = 0. For fi > /x c , the color symmetric phase becomes 
unstable because a solution with A = does not minimize f2. Here, the solution with to 7^ and A ^ 0, corresponding 
to the 2SC phase, gives the global minimum of £1, and thereby the color superconducting phase is favored (note that 
in the 2SC phase fi > to, whereas /i < to in the color symmetric one). The transition between these two phases is of 
the first-order, which is characterized by a discontinuity in the behavior of m vs. fi (see Fig. 1). 



We are now interested in the investigation of the modification of meson and diquark masses in dense and cold quark 
matter with the color symmetry broken because of the 2SC diquark condensation. In the vacuum (T = 0, fx = 0), 
particle masses are obtained from propagator poles, or alternatively, from zeros of the one-particle irreducible (1PI) 



An expression for the thermodynamic potential for a system of free fermions can be found elsewhere; see, e. g., 




(17) 



(18) 



III. MESONS AND DIQUARKS IN DENSE QUARK MATTER 
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two-point Green's functions. Since the Lorentz invariance is preserved in the vacuum, these functions in (Minkowski) 
momentum space depend on p 2 = p\ — p 2 only. The squared mass of the particle is thus equal to the value of p 2 
where the corresponding two-point 1PI correlators (Green's function) vanishes. Insofar as p 2 is a Lorentz invariant, 
one can, for simplicity, choose the rest frame (p = 0), put p 2 — p^ and then consider these correlators as functions of 
Po alone. On the contrary, in a dense medium, the Lorentz invariance is broken, so the two-point Green's functions in 
momentum space should be treated as functions of two variables: po and p 2 (the rotational symmetry is presumably 
conserved). The zeros of 1PI two-point Green's functions in the po plane will determine the particle and antiparticle 
dispersion laws, i. e. the relations between their energy and three-momenta. In this case, the scalar particle mass is 
defined as the value of the particle energy at p = (see, e. g. 25]). Recently, the Bethe-Salpeter equation approach 
has been used to obtain diquark masses in the 2SC phase of cold dense QCD at asymptotically large values of the 
chemical potential 26]. There, the mass of the diquark was defined as the energy of a bound state of two virtual 
quarks in the center of mass frame, i. e. in the rest frame for the whole diquark. Let us note here that just this 
quantity is measured on the lattice, where it is given by the exponential fall-off of the particle propagator at large 
Euclidean time (see, e. g., [23] )- 8 Similarly as in [2^, |2^] and in numerous other papers, we shall denote, in the 
following, the rest frame energy of a composite scalar particle (meson or diquark), moving in a dense medium, as 
"mass". (In general, the values of the mass depend on the chemical potential.) 

Any 1PI Green's function can be found from the effective action S e g, which up to second order in boson fields has 
the form 

S ° s = Ijl I d 4 xd 4 yX(x)Tl XY (x,y)Y(y) + • • • , 

1 X,Y 

where X, Y — ir a , er, A^, A^, and Hxy( x , y) is the coordinate representation of the 1PI Green's function for the fields 
X, Y . Instead of employing functional integration to get S c ([ and then Il^r, we shall use implicitly Feynman diagrams. 
Starting from the Lagrangian ©, one can expand the resulting effective action S e s in a power series of meson and 
diquark fluctuations. Keeping there only the second order contributions, we immediately obtain the 1PI two-point 
correlators in the one-loop approximation. By considering this functions as functions of po at zero three-momentum 
(p = 0), we shall next analyze their zeros which, as explained above, will give the masses of resonances. 



A. Pion mass 



Let us begin with the calculation of the pion mass. In the momentum representation, the 1PI two-point function 
II T r a ^ i) (P) for the pion has the following form (all calculations are performed in the rest frame, P = (po, 0, 0, 0) ) 

IWP) = + \j ^Tr [S{P + q)TlS{q)Ti] . (19) 

Here, the vertex of the pion-quark interaction is given by the 2x2 matrix = ( *^ 5Ta 9 t ) . The first term in 



, i^Tl 

the right hand side of l|19|) is the tree contribution from Lagrangian 10, the second term arises from the one- loop 
diagram with two pion legs. After intermediate trace calculations, this function takes the form 

Sab , f d 4 q qoipo + Qo) -E+E- - A 2 



n w . W6 (P) = -^- + i6t<5 o6 J 



(2tt) 4 D-{qo)D + (p + qo) 



^s-A f d ' q go(po + go) -E+E- 

+mab J (2tt)4 [(p + qo )2 - (E+n(q ) 2 - ( E -) 2 ] " ^ 

Implementing the imaginary-time formalism, as described in the end of the previous section (see also Appendix A for 
details), we reduce eq. ()20fl to three-dimensional integration in momentum space 



iWPH-ll + s.u J 



Sab , „ f d 3 q E+E A + E+E- + A 2 E++E, 



A 



(2tt)3 E+E a (E+ + E A ) 2 -p 2 



'Wf%^^ S ^M. (21) 



In the last case the lattice calculations were performed in some QCD-like theories (QCD with two colors, etc), where the fermion 
determinant is positive even at nonzero /i. Evidently, this mass definition is borrowed from particle physics. Note that in condensed 
matter physics usually the term "energy gap" is used for this quantity (see, e. g., L'lJ . where both terms, "mass" and "energy gap", are 
used for the rest frame energy of scalar particles). 
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FIG. 2: The masses of the cr-meson (solid line) and pion (dashed line) as functions of ji. 



Then, up to a sign, the pion unnormalized propagator equals (n 7ra7rb (P)) _1 , whose pole in p$ is given by the zero of 
the function Il W7r (po) from l|21|l . We have searched for the roots of the equation II ffW (po) = numerically; the results 
for the pion mass M n at various fj, are plotted in Fig. 2. A similar behavior of the pion mass in the Cooper pairing 
phase of a dense quark matter was found in [2^| in the framework of a two-colored NJL model. 

In the color symmetric phase, the pion mass is below the threshold for the pion decay to a quark-antiquark pair, 
and the pion, therefore, is almost stable (only electroweak decay channels are allowed). Moreover, we have found that 
in the 2SC phase the pion is also an almost stable particle. This conclusion is supported by the following arguments. 
On the one hand, the pion cannot decay into a red/green quark-antiquark pair, since M T is less than the minimal 
energy E m i n = \J1jjl — to) 2 + |A| 2 + \f(p + to) 2 + |A| 2 that is necessary to create these particles (evidently, E m [ n is 
the value of E^ + E^ (see eq. (|12fl and the text after it) at E = to). On the other hand, the decay of a pion into a pair 
of blue quark-antiquark is also forbidden: i) Since all states with E < /x are occupied by quarks from the medium, the 
pion cannot decay into a blue quark with such values of E because of the Pauli blocking, i.e. the analogy to the Mott 
effect, ii) In the region E > //, the energy of a blue quark-antiquark pair (E~ + E + ) takes its least value -Emm = 2/i 
at E = /i, so the pion decay is again impossible, since < 2fi. As a consequence, the pion is almost stable in the 
2SC phase, too. 



B. Mixing of 0--A2 in the 2SC phase. Scalar meson mass 

An investigation of the 1PI Green's functions H a x(P) with P = (po, 0,0,0) for the fields X = A* A , As shows that 
in the 2SC phase the cr-meson is mixed with the A2 diquark. (At /i > /i c , such a mixing occurs in the NJL model with 
two-colored quarks 28], too. Moreover, as our preliminary results show, the mixing is present even if the condition 
of color neutrality of the 2SC phase is imposed.) The mass of the cr-meson in this case is given by the solution of the 
equation det(II) = where II(P) is the 3x3 matrix 

(n aa {p) , iwp), n CTA| (P) \ 
n(p)= n A2ff (p), n A2A2 (p), n AaA .(p) . (22) 
Vn A « CT (p), n A - A2 (p), n A , A .(P)/ 
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(Up to a sign, II(P) is the inverse propagator matrix for the cr-meson and A?j, A2 diquarks.) After tedious but 
straightforward calculations, similar to those in the pion case, we get 

1 , 1fiA 2 2 f A f 1 1 



(2 7 r)^+[ p 2_ 4(jE j+ ) 2 ] ' •■ 7 (27r)3£ A [ p 2-4(£ A ) 



(27r)3^ ^+[4(^)2-^] £ ;- [ 4(£-)2_ p 2 

/• ri 3 g g E+E A A 2 ^ + ^a . 1fi f d 3 q <F 9{E-n) 
J (2^ E* E + E - (E++E A r~ P y J (2tt)3 £ 4£ 2 -p 2 ' I ; 

n CTA2 (p) = n A . ff (P) = n CTA .(-p) = iw(-p) - 

= 4toa/^L/ 2E++P0 + 2E ~- P0 1 (24) 

n A2 A 2 (P) - n A|A .(P) = 4A 2 i q (pD, n AaA .(P) = ii a » A2 (-p) = 

= - i ^ + / A + (4A 2 -2p 2 )/ (^) + 4 p / 1 (^ )j (25) 

where J A is given by 0, Io(Po) = ^+ + ^-j h(jPo) = B + - B_, and 

/■ d 3 g 1 , f d 3 q 1 ^ 

S+= 7 (2JrE+[pl-A(E+rr B - = J (2nfE A [pl-4(E A n (2?) 

(The formulas (|23fl -H25 [) are valid both in the 2SC (A ^ 0) and in the color symmetric (A = 0) phases. In the 2SC 
phase, the 1PI Green's functions II A2A j(P), II A * A2 (P) (|2"5)l get a more simple form if one uses the identity 1 = 4G2IA 
following from the gap equations (see eq. |T%)l).) In the general case (m / 0, A ^ 0), the equation det(II) = has 
a rather complicated form. Fortunately, in the color symmetric phase with A = 0, i. e. at fi < fi c , the nondiagonal 
terms Ii a x{P) (X — Ajjj, A2) in ^22(1. which are responsible for the mixing of the cr-meson and the diquark, vanish 
because they are proportional to A. Therefore, the cr-meson mass decouples from the diquark spectrum and is found 
from the equation 

IWP) = 0. (28) 

On the other hand, in the 2SC phase (see Fig. 1), the constituent quark mass to is small (or even equal to zero if 
too = 0) in the 2SC phase (see Fig. 1), so one can ignore the nondiagonal elements II (tA2 (P), n CTA *(P) in II because 
they are negligibly small (as proved a posteriori by numerical computations), and the cr-meson mass M a is again 
found from eq. (|28H . The numerical solution of eq. I|28|) is presented in Fig. 2, where one can see that both a- and 
7r-meson masses are increasing functions of /1 in the 2SC phase. 9 At the same time, the difference between M a and 
M„ decreases with p,; 5M = M a — M w becomes negligible at sufficiently high /i, which is understood as an evidence 
of the chiral symmetry restoration. The decrease of the dynamical quark mass to at large /j, (see Fig. 1) is also in 
accordance with this conclusion. 

Finally, we would like to note that the cr-meson is almost stable in the 2SC phase for the same reason which was 
explained in the last paragraph of the previous section for the case of the pion. 

C. Diquark masses 

In dense quark matter (at nonzero (baryon) chemical potential), the symmetry of the Lagrangian under charge 
conjugation is violated by the chemical potential. As a consequence, the mass spectrum of diquarks can split, and 
diquarks will differ from antidiquarks not only by charge but also by mass. 

1. Diquark masses in the color symmetric phase (A — 0) 

In the color symmetric phase (/i < fi c = 350 MeV), the ground state of the quark matter is described by A = 0, and 
there is no mixing between diquarks in the one-loop approximation. Therefore, one can, e. <?., consider the propagator 



9 Note that the term in 1231 that is proportional to A 2 m 2 is comparable or even less than nondiagonal elements H a x(P) (X = A|, A2). 
So it was neglected in numerical analysis of 1281 . 
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of A£, A2 alone. It follows from eq. (|25(l that Ha 2 a 2 (P) = Ha*a*(-P) = at A = 0, and one needs only 

iWP) = n A3 A 2 (-P) = - 4^ - 16 / (0^t^i^ - -4^ + (M) 

Here, P = (pg, 0,0,0), e — (pq + 2/i) 2 . (In obtaining l)29|) . we have used the relation /j, < m, i. e. E~ > 0, which is 
true for the color symmetric phase.) Then, the 2x2 inverse propagator matrix G^iP) m the A2,A 2 -sector of the 
NJL model has the form 

s "< p '=-(n 4 ; 4 >), n4 ' 4 o s(p) )' (30) 

Clearly, the mass spectrum is determined by the equation det(^ 2 ~ 1 ) = IIa 2 a* (-P)IIa* A 2 (-P) = 0, or by zeros of l|29|l . 
where the function F(e) is analytical in the whole complex e-plane, except for the cut 4m 2 < e along the real axis. (In 
general, the function F(e) is defined on a complex Riemann surface which is described by several sheets. However, 
a direct numerical computation based on eq. Q29[) gives its values in the first sheet only. To find a value on the rest 
of the Riemann surface, a special procedure of continuation is needed.) The numerical analysis of (|29|l in the first 
Riemann sheet shows that the equation IIa 2 a*(-P ) = Ha*,a 2 (— P) — has a root (eo) on the real axis (0 < eo < Am 2 ), 
providing us with the following massive diquark modes 

M A = 1.998m - 2/i, Ma* = 1.998m + 2/i. (31) 

We relate Ma in l|31|) to the mass of the diquark with the baryon number B — 2/3 and Ma* to the mass of 
the antidiquark with B — —2/3. (Qualitatively, a similar behavior of diquark and antidiquark masses vs. /< was 
obtained in [2^ in the NJL model for two-colored quarks.) It follows from (|31|l that in the vacuum (/t = 0) the 
diquark/antidiquark mass is ~ 2m. Clearly, in the color symmetric phase at /i ^ both quantities Ma and Ma* 
from 1131 [I are nothing else but the rest frame excitation energies for a diquark and antidiquark, respectively. 

In the color symmetric phase for /t < m the diquarks and antidiquarks are stable. Indeed, the diquark mass Ma is 
smaller here than the energy 2(m — \x) necessary to create a pair of free quarks. Finally, due to the underling color 
SU(3) C symmetry, the previous statement is valid also for Ag, A5 and A£, A7. As a result, we have a color antitriplet 
of diquarks with the mass Ma (|31|l as well as a color triplet of antidiquarks with the mass Ma* . The results of 
numerical computations are presented in Fig. 3: the solid line shows the behavior of the antidiquark triplet mass Ma* 
in the region of fi < /i c — 350 MeV whereas the dashed line corresponds to the antitriplet. 

In our analysis, we have used the constraint G 2 = 3Gi/4, thereby fixing the constant G2 through G\. It is useful, 
however, to discuss now the influence of G2 on diquark masses. Indeed, it is clear from 129|) that the root eo lies inside 
the interval < e < 4m 2 only if G 2 < G 2 < G 2 *, where G 2 and G 2 * are defined by 



Gl EE 
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4F(4m 2 ) 4 [A %/m 2 + A 2 + m 2 ln((A + Vm 2 + A 2 )/m)] ' 



G? = — = - = 3mGl ■ (32) 

4F(0) 4 [AVm 2 + A 2 - m 2 ln((A + Vm 2 + A 2 )/m)] 2(m - m ) 

In this case, there are stable diquarks and antidiquarks in the color symmetric phase. The behavior of their masses 
qualitatively resembles that given by eqs. 1|31|) . For a rather weak interaction in the diquark channel (G 2 < GJ), eo 
runs onto the second Riemann sheet, and unstable diquark modes (resonances) appear. Unlike this, a sufficiently 
strong interaction in the diquark channel (G 2 > G 2 *) pushes eo towards the negative semi-axis, i. e. (po + 2/i) 2 < 0. 
The latter indicates a tachyon singularity in the diquark propagator, evidencing that the color symmetric ground state 
is not stable. Indeed, as it has been shown in |30j . at a very large G 2 the color symmetry is spontaneously broken 
even at a vanishing chemical potential. 



2. Diquark masses in the 2SC phase (A/Oj 

Let us now focus on the masses of A^, A 2 -fields. As we have already shown in Section Till Bl these diquarks are 
mixed with the cr-meson in the 2SC phase because of the nonvanishing terms ilg^ (P) an d n^Aj (P) m the matrix 
n(P) 122(1 . However, keeping in mind that the constituent quark mass m is small in the color superconducting phase, 
one can ignore this mixing. The problem becomes, thereby, drastically simplified, and one just has to calculate the 
determinant of 



a-u P \ -_( n A 2 A 2 (n n A2 A*(P) \ 
y " { >~ l,n A . Aa (p), n A . AS (P) J 



(33) 
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FIG. 3: The masses of diquarks. At fi < [i c — 350 MeV, six diquark states are splitted into a (color)triplet of heavy states with 
mass Ma* (solid line) and an anti-triplet (dashed line) of light states with mass Ma- In the 2SC phase (n > fj, c ), one observes 
3 massless diquarks (dashed line): a doublet of light diquarks with the mass Mn g ht (dotted line) and a heavy singlet state with 
the mass M (dash-dotted line). The shaded rectangular displays the width of the heavy singlet resonance; its upper border is 
half-width higher than the mass and the bottom border is half-width lower. 



and equate it to zero. The resulting equation will determine the masses of Ajjj, A2. Taking the expressions for the 
matrix elements in l|33[) from eqs. (|24|l and i|25|) and using the relation 1 = (valid in the 2SC phase only), we 

get the mass equation 

det^ 1 ) = Ap 2 {{p 2 - AA 2 )I 2 (p 2 ) - 4Z?(p§)} = 0. (34) 

It has the apparent solution p\ — corresponding to a Nambu-Goldstone (NG) boson. 10 The second solution of 1341) 
exists on the second Riemann sheet for p\ only (see Appendix B). 
Near a zero, the determinant det(5^~ 1 ) can be approximated by 

det(g 2 _1 ) ~ p\ - M 2 + iMT. (35) 

Here, M is the mass of the resonance, and T is its width. Let po be a root of the equation 

Pi - M 2 + iMT = 0, (36) 

the mass and width are then given by 

m = v /r^, r = -^. (37) 

For a small width, we can write for the root of eq. (|36() 

ffcwM-i-. (38) 

The appearance of the imaginary part in (|38() is a consequence of the fact that this diquark is a resonance in the 
2SC phase and can decay into free quarks. Our numerical estimates for M and T at various \i are plotted in Fig. 3. 



The equation det(II(po)) = where II(po) is given by 1221 also has a NG-solution (p^ = 0) in the 2SC phase. Indeed, one can easily see 
that the elements of the second and third columns of II(po) are equal at po = 0, and the determinant of this matrix is thereby equal to 
zero at po = 0. Since det(II(po)) is an even function of po, one can conclude that this solution is doubly degenerated. 
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The dot-dashed line corresponds to M, while the width is given by the hight of the shaded block, half-width up and 
half-width down. 

We would like to point out here that both the resonance and the above mentioned NG-boson are color singlets 
with respect to SU(2) C . Since the obtained mass M is much greater than even twice the energy y/(m + fi) 2 + \ A\ 2 
necessary for creating an antiquark in the 2SC phase, there is no wonder that this diquark mode is unstable, unlike 
the pion and tr, which are stable due to the Mott effect. 

A detailed investigation of the diquark masses in the AgjAs- and Ay,A7 sectors has already been done in our 
previuos paper |jj . It was found there that in each of these sectors there was an NG-boson as well as a light diquark 
excitation with the same mass Mi; g ht proportional to the color-8 charge density of the quark matter in the 2SC 
ground state. So, we can conclude that, in total, there is an abnormal number of three NG-bosons in the theory, 
instead of expected five (Since SU(3) C is spontaneously broken down to SU(2) C , five group generators affect the 
diquark condensate, and therefore, five NG bosons should appear, according to the Goldstone theorem). This entails 
absence of baryon superfluidity in the 2SC phase (see also the discussions in [2!j[j^]], where the effect of NG-boson 
deficiency was observed in other relativistic models with broken Lorentz invariance). The light diquarks form a stable 
SU(2) c -doublct with the mass Mu g ht whose dependence on [i is also shown in Fig. 3 (dotted line). 

IV. SUMMARY AND DISCUSSION 

In the present paper, we have investigated the mass spectrum of meson and diquark excitations in cold dense 
quark matter. We started from a low-energy effective model of the Nambu-Jona-Lasinio type for quarks of two 
flavors including a single quark chemical potential, for simplicity. Despite of the lack of confinement, this model quite 
satisfactorily describes the masses and dynamics of light mesons in normal quark matter (at rather small values of 
chemical potential). Since the investigation of color superconductivity became popular nowadays, NJL models have 
also been widely used to explore, in particular, the quark matter phase diagram for intermediate densities, i. e. under 
conditions where all other approaches fail. 

Using the one-loop approximation, we calculated two-point correlators of mesons and showed that the masses of 
7r- and a- mesons grow with the quark chemical potential in the 2SC phase (see Fig. 2). The mass difference between 
them vanishes at asymptotically large fi, in accordance with chiral symmetry restoration. Moreover, these mesons 
are stable in the 2SC phase due to the Mott effect. As far as we know, the properties of it- and a-mesons in the 2SC 
phase have not been discussed in the literature before. 

In the diquark sector, the situation is more involved in the 2SC phase. Indeed, when the color SU(3) c -symmetry 
is spontaneously broken down to SU(2) C , one naturally expects five (massless) NG bosons to appear. Unlike this, 
one finds only three massless bosonic excitations Q: a color singlet and a color doublet (due to the residual SU(2) C 
symmetry) . In spite of the abnormal number of NG-bosons (notice that each member of the doublet has a quadratic 
dependence of its energy on three-momentum when it is almost at rest), there is no contradiction with the NG-bosons 
counting [32|. Apart from this, there are also two light and one heavy diquark modes (see Fig. 3). The first two 
are stable, whereas the last one is a resonance with finite width, and its 1PI Green's function possesses a zero in the 
second Riemann sheet for the energy variable. 

We have also found that the antidiquark masses exceed those of the diquarks in the color symmetric phase (for 
/i < /j, c = 350 MeV). This splitting of the masses is explained by the violation of C-parity (charge conjugation) in the 
presence of a chemical potential. In contrast, at /i = the model is C-invariant and all diquarks and antidiquarks 
have the same mass which is slightly lower than two dynamical quark masses, ~700 MeV. Our result for the diquark 
mass in the vacuum (/x = 0) is in agreement with |33J, where the value as large as ~ 800 MeV was claimed to follow 
from QCD via solving a Bethe-Salpeter equation in the rainbow-ladder approximation. 

Of course, all observable particles render themselves as colorless objects in the hadronic phase, and the diquarks 
are expected to be confined, as they are not SU(3) C color singlets. Nevertheless, one may look at our and other 
related results on diquark masses as an indication of the existence of rather strong quark-quark correlations inside 
baryons, which might help to explain baryon dynamics. Some lattice simulations reveal strong attraction in the 
diquark channel [34j with a diquark mass ~600 MeV. Recently, in [35J, the mass and extremely narrow width, as well 
as other properties, of the pentaquark + were explained just on the assumption that it is composed of an antiquark 
and two highly correlated ud-pairs. At present time, the nature of the mechanism which may entail strong attraction 
of quarks in diquark channels is actively discussed both in the nonperturbative QCD and in other models (see, e. g. 
|3("i| and references therein). 

Finally, let us comment on the fact that we considered in this paper only a single chemical potential /i (common for 
all quarks). Obviously, in such a simplified approach, the 2SC quark matter is neither color nor electrically neutral, 
as it would be expected for realistic situations like, for example, in the cores of compact stars. In order to study color 
superconductivity in the case of neutral matter, one has to study more complex NJL models including several new 
chemical potentials 5]). Despite this drawback, the chosen simplified approach, nevertheless, seems to us interesting 
enough to get some deeper understanding of the dynamics of mesons and diquarks in the color symmetric and 2SC 
phases of cold dense quark matter. A generalization of this approach to more sophisticated NJL models including 
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several chemical potentials is presently under investigation. 
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APPENDIX A: LOOP INTEGRALS ATT/0 AND fj, ^ 

To evaluate loop integrals, we use in our paper the imaginary-time formalism (see, e. g. p^). First of all, it is 
supposed that the system is in the thermodynamic equilibrium with a thermal bath of some temperature T. As 
usual, to study a hot and dense system, one has to calculate thermal Green's functions (TGF), which are periodic 
(for each boson field) and antiperiodic (for each fermion field) functions of imaginary time. Their Fourier images 
are, thereby, functions defined at discrete points in the imaginary axis po = iuj n , n — 0, ±1,±2, . . ., with uj n being 
Matsubara frequencies. When calculating a loop diagram, one has to sum over qo = iuj n = i2mrT for bosons and 
<7o = i^n = i(2n + l)irT for fermions. In accordance with this, the integration over go in (f 141) . (|15fl . etc is to be 
performed in two steps: First, all momenta are considered in Euclidean metric, and a sum over Matsubara frequencies 
is performed. Then, the limit T — * is reached and the resulting TGF are continued in the complex plane to the real 
axis, which corresponds to their definition in Minkowski metric. As a result, we have then for 1PI Green's functions 
the formulas where the four-dimensional integration is reduced to three dimensions. 

Let us explain the scheme, briefly outlined above, in a bit more details. First, we delay the three-dimensional 
integration until the end of the calculations, and consider only one-dimensional go-integrals. Only one such integration 
comes from each diagram in the one-loop approximation for the 1PI Green's functions with two bosonic legs, and it 
can be represented as 

F(po)= I ^/(MPo), (Al) 
J 2tti 

where the function f(qo-,q',Po) is a product of vertices and quark propagators from a one-loop diagram, while the 
external three-momentum in each diagram is supposed to be zero, p = 0. The dependence of the function F(po) on 
the momentum q is implied, however, we do not write q explicitly, since only the dependence on po is important for 
the moment. 

At T / 0, one would obtain in the right hand side of (|A1() a sum over the fermionic Matsubara frequencies 
uj n = (2n + l)wT instead of the integral. Let us denote the result of the sum by the function Ft(pq). It is defined 
only at the values of po — iv^ = i2knT, corresponding to bosonic Matsubara frequencies, because each external line 
in the implied one-loop diagram corresponds to a boson. The result looks like 

oo 

F T (iv k )=T ^2 f(iuj n ,q;iv k ). (A2) 

n— — oo 

Let us assume that the function / (u> , q; iv k ) falls sufficiently quickly in any direction on the complex w-plane and 
contains only simple poles everywhere, except for the imaginary axis. In this case, using the Cauchy theorem, one 
can replace the sum over Matsubara frequencies in (|A2|) by " closed" contour integrals and obtains 

oo If duj CO 

F T {iv k ) = T^2 f{iu n ,q;ivk) = - j> —f(u,q\iv k )taxih^—^, (A3) 

n ="°° c 1+ c 2 

where the contour C\ is just the straight line from — ioo + e to ioo + e, and Ci is the straight line from ioo — e to 
— ioo — e. The contour integral in l|A3(l is indeed the sum of two integrals: C\ and C^- Both C\ and Ci can be closed 
by infinite arcs in the right and left halves of the complex plane because the integration along these arcs vanishes if 
the integrands fall quickly enough near infinity. As a consequence, we can then integrate along the closed contours 
C\ and Ci. As the integrand was supposed to have only simple poles, the Cauchy theorem immediately gives us 
the result of integrations via a sum of residues of the integrand at these poles, which are determined by the poles of 
the quark propagator. It is important to keep in mind before taking the limit T — > that v k — 2knT. The point 
is that after the calculation of all residues the frequency iv k will appear in the argument of the hyperbolic tangent 
tanh(w/2T), which is periodic on the imaginary axis and the period is just equal to i2?tT. Therefore, the tangent does 
not depend on iv k , and one can put i>k = in its argument. After this, one can reach the limit T — > and continue 
the function Ft\t=o to real energies p$, which is formally obtained through the substitution v k — > —ipo in the rest of 
the expression. 
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We follow the scheme explained above to calculate the tadpole contributions (see l|13|) ). which do not depend on 
external momenta. Let us consider, as an example, the one-loop tadpole contribution for (A^x)), which is nothing 
else than the right hand side of l(TK|) . Replacing the qo-integral by a sum over Matsubara frequencies w n = (2n+l)nT 
(see UlSf) 1 ) and following other steps, we obtain 



1 



D-(icj n ) 



dui , / uj \ f d 3 q ( 1 1 



= -2 t —^-)J—\— + > , (A4) 



Ci+c 2 



where D±(ui) arc defined by l|12l) . and the two clockwise contours C\ and Ci enclose the right and left halves of the 
complex w-plane. The integrand in l|A4|l has only four simple poles at uj = and ui = ±E^. Finally, we replace 

the integral in (|A4|) by a sum of the residues at these poles, multiplied by 2iri (according to the Cauchy theorem), 
and take the limit T — > 0. This gives for I a = limy^o Ia{T) the expression displayed in the right hand side of (|TH|) . 



APPENDIX B: SEARCHING FOR THE RESONANCE SOLUTION OF EQ. iHTHl 

The nontrivial solution of Ij34|l obeys the equation 

(z/ -2/i)(z/ + 2/i) = 0, (Bl) 

where z 2 = p\ — 4A 2 and 7 .i are given by (|25|) . Below we shall ignore, for simplicity, the dynamical quark mass 
(i. e. putting m = 0), because we assume (and this assumption is a posteriori corroborated by numerical calculations 
of the heavy diquark mass both for m = and m 7^ 0) that this simplification does not strongly affect our results. 
First of all, we integrate over the angles in (|26|) and l|27() and introduce a new variable y = (E + fi) 2 , instead of the 
three-momentum, in the integrals A+ and -B+ (recall that E = \q\) 

A+ = J_ (Vv-rfdy B+ = — ( (A+ " ? WlLjWv (B2) 



Since 



A_ = f \.-.)dE+ [ {■■■)dE, (B3) 

JO J [i 

one can introduce the new variables yfy = \i — E and y/y = E— fi for the first and second integrals of (|B3|I . respectively, 
and get 

47r 2 i ^/yy / y + A 2 [z 2 -Ay] 4tt 2 J y /y^/y~+~K 2 [z 2 - Ay] 

and 



Jo _ + A _ _ 1 (VV-^v + 1 (Vv + ,)^ 

^ Jo Jyy/y + A 2 [z 2 - Ay] ^ Jo Jyy?y + tS. 2 [z 2 - Ay] 



In a similar way, one obtains 



, R R 1 /- (A+A1)2 (Vy-v) 2 dy 1 /- (A ^' )2 (Vy + v) 2 dy 



4tt 2 7 y / y + A 2 [z 2 -4y] 4tt 2 7 V7+A 2 [z 2 - 4y] 

The quantities /o,i are analytical functions of the complex variable z 2 = a — ib on the whole complex plane, except 
for the cut L on the real axis, defined by < z 2 < (A + pi) 2 (the first Riemann sheet). A numerical processing of 
the integrals (|B5|I and (|B6(1 gives the values of functions Iq and I\ on the first Riemann sheet only. But there is no 
solution for eq. (|B1|I in the first Riemann sheet because the root of eq. IjBlfl lies on the lower half-plane (b > 0) of the 
second Riemann sheet, and in order to find it, we have to continue Iq and I\ to the second sheet. When z 2 crosses 
the cut L downwards from the upper half of the first sheet, the integrals Iq (JB51 and I\ (|B6|I become singular because 
of the apparent poles in the integrands. Since an integral does not change when the contour is being continuously 
transformed on the complex plane until it crosses a singularity of the integrand, we carefully deviate our contour on 
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(A±m) 2 



C'z 




FIG. 4: The deformation of the integration contour on the complex j/-plane for the functions when the parameter z 2 

moves onto the lower half-plane of the second Riemann sheet. 



the complex y-plane to circumvent the singularity in a way shown in Fig. When the contours C\ and C2 overlap 
each other, the sum of the integrals along them vanishes, and as a result, we obtain that Iq and I\ in the lower 
half-plane of the second Riemann sheet differ from i|B5|) and l|B6|) by an additional integral over the contour C around 
the pole z 2 (see Fig.0J which is equal to the residue of the integrand at z 1 . Let us denote the continuations of Iq and 
Ii to the second Riemann sheet as Iq and I\. Then, one has 

i(z 2 + 4u 2 ) ~ izn 

I =I V ; , h=h + ; (B7) 

iirzVz 2 + 4A 2 2Wz 2 + 4A 2 V ; 

where z = \[~z 1 ~ = a\ — ibi, with both a% and b\ being real and positive: a\ = \l + a 2 + S 2 ", &i = 
bj\j2a + 2\/a 2 + fr 2 ; apart from, V^ 2 + 4A 2 = a x - ibi, where fii = ^/±(a + 4A 2 ) + ±i/(a + 4A 2 ) 2 &i = 

Vy / 2(a + 4A 2 ) +2v / (^T4A 2 ) 2 + fe 2 . 

Let us now consider the function F{z 2 ) = \fz? t, Iq{z 2 ) + 21\{z 2 ) (it is the second multiplier in IjBTfl ). Its continuation 
to the second Riemann sheet is F — > F{z 2 ) — \/z^I a + 2I\ = F(z 2 ) — {(Vz 2 — 2/i) 2 /(47r\/z 2_ T4A 2 ). Numerical 
solution of the equation F(z 2 ) = at various fi > fi c gives one root per one value of /i: z 2 + 4A 2 = M 2 — iMT. For 
example, if one puts A = 115 MeV, fi = 350 MeV and A = 618 MeV, one gets the mass M = 1111 MeV and the 
width T = 446 MeV. 
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